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Abstract 

CERN LHC provides a good experimental platform to perturbatively probe the fundamental 
gravity scale up to several TeV, with the precise value depending on the number of extra dimen- 
sions. The leading experimental signal of graviton at LHC is from the process pp — > jet + flj,, where 
fij, stands for the transverse missing energy. A detailed discussion on the hadronic production of 
real graviton through hard subprocesses: qq^G + g,g + q^G + q and g + g — > G + g have been 
studied within the quantum gravity theory with large extra dimensions. The main theoretical 
uncertainties together with the dominant standard model background to these processes, e.g. 
qq — > Z° + g and g + q — > Z° + q with Z° further decaying into neutrinos, have also been discussed. 
It is found that only in certain jet energy region and with certain number of extra dimensions 
can the quantum gravity signal be distinguished from the background, which inversely lead to the 
effective scale to be probed up to (8.8 ± 0.9) TeV for two extra dimensions, and (5.9 ± 0.5) 
TeV for four extra dimensions with sufficient integrated luminosity, e.g. 100/6 -1 , at CERN LHC. 
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I. INTRODUCTION 



It has been purposed that standard model (SM) particles live in the usual 3 + 1- 

n n 

dimensional space, while gravity can propagate in a higher-dimension space PJ, 12|]. Such 
a scenario is helpful to reduce the fundamental mass scale from the large Plank scale down 
to be about TeV-sc&le and then to solve the so-called hierarchy problem, so it arouses peo- 
ple's interests since its first announcement. Numerous attempts have been carried out to 
find the signal of the large extra dimension, i.e. measurements of gravity at short distances, 
studies of various astrophysical and cosmological implications of large extra dimension and 



n n 



collider searches for virtual and real graviton effects [3j, llO|]. The CERN 

LHC, with its high collision energy (14 TeV) and high luminosity (10 34 cm _2 s _1 ), shall pro- 
vide a better platform to study the extra-dimension phenomenology both experimentally 
and theoretically. 

The leading experimental signal of real graviton at LHC is from the hadronic process 
PP ~ > J e ^ + fir w hh ftp stands for the transverse missing energy. We shall present a detailed 
discussion on the hadronic production of graviton through the hard subprocesses: qq — > 



G + g, g + q — > G + q and g 



large extra dimensions 
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g — > G + g. By using the quantum gravity effective theory with 



121 ]. we make a try to study the quantum gravity effects and its 



dependence on the number of extra dimension, and to study up to what energy scale can 
LHC probe. Furthermore, main theoretical uncertainty for the graviton production shall be 
studied, which includes the parton distribution functions (PDFs) for the initial partons, the 
choice of the factorization scale pL 2 F and the typical energy scale Q 2 for the hard scattering 
amplitude, the number of the extra dimension for the graviton production, and etc.. 

The processes qq —>■ g + Z° and qg — > qZ°, followed by an invisible decay of Z°, i.e. 
Z° — > vv y give an irreducible physical background to graviton production. We refer to 
these processes as the 'SM background'. We will estimate the observability of the graviton 
signal by comparing its hadronic cross-section to that of 'SM background'. There are other 
important background sources from miss-measured jets and W production with forward 
leptons, however as argued in Ref . [3| that these backgrounds decrease sharply as the lower 
bound on missing Et is increased, so we shall not considered it here. 

The remainder of the paper is organized as follows. Section II is devoted to give the main 
formulae for graviton production within the framework of the quantum gravity effective 
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theory with large extra dimensions. Numerical results and discussions are presented in 
Section III, where the uncertainties in estimates and a discussion on the value of the effective 
energy scale Mp> shall be presented. The last section is reserved for a summary. 



II. CALCULATION TECHNOLOGY 

According to the QCD factorization formulae, the hadronic production of graviton at the 
collision center of mass energy can generally be written as 

da(S,E T ,---) = 

12 f f dx 1 dx 2 F} Il Pl (xi,ii%)Fj l2! p 2 (x2,iJ^)dai j ^j et ^^ ,(1) 

where F % Hi Pi (xi,fi F ) and F 3 H p 2 (x 2 , (i 2 F ) are PDFs of incoming hadrons Hi (momentum P x ) 
and if 2 (momentum P 2 ) for parton i (with the momentum fraction x\) and parton j (with the 
momentum fraction x 2 ) respectively. Q 2 is the 'characteristic energy scale of the subprocess 
squared'; and /i F is the 'energy scale squared' where the factorization about the PDFs and 
the hard subprocess is made. Usually for leading order (LO) calculation to obtain the best 
results, the two scales \xp and Q 2 are carried out as the same, thus later on we take fi F = Q 2 
except one case when estimating the uncertainty from LO and the ambiguity of the choices 
about fi 2 F and Q 2 . 

d&ij^jet+p^ stands for the differential cross-section of the relevant hard subprocess, in 
which s = X\X 2 S is the center of mass energy of the subprocess and stands for the 
missing transverse energy. Within the framework of the quantum gravity theory with large 
extra dimensions, the differential cross-section for inclusive graviton (G) production, i.e. 
ij — > Gk with i, j and k stands for corresponding partons, can be written as [3] 

d 2 « Q Mp J _ 1 da m 



dtdm l3 ^ k ~ M 2 D +S dt 



where da m /dt stands for the differential cross-section for producing a single Kaluza-Klein 
graviton of mass m, M P = M p /\/8tt = 2.4 x i0 18 Ge^, Ss-i is the surface of a unit-radius 
sphere in 5 dimensions, for 5 = 2n and n integer, S$-i = 2-n n jin — 1)! and for 8 = 2n + 1, 
Ss-i = 2tx ti /Yi^Z^ik + |). For the present 2 — > 2 subprocesses, the Mandelstam variables are 
defined as s = (p« + Pj) 2 , t — (pi — pc) 2 and u = (pi — Pk) 2 - Since the graviton interaction 
vertex is suppressed by 1/M P , it can be found that a m oc Mp 2 , and the factor Mj, appearing 



from the phase-space summation exactly cancels the Planck mass dependence in Eq. (J2j). 
In another words, the large phase space of the Kaluza-Klein modes, corresponding to the 
large volume of the compactified space, exactly cancels the dependence on Mp and gives an 
effective interaction suppressed only by inverse powers of Mp>- This is the reason why sizable 
contributions from the graviton may be observed at LHC, and inversely, one can estimate to 
what energy scale can LHC probe. For the differential partonic cross-sections da m /dt that 
produce a single Kaluza-Klein graviton with mass m, we obtain 

^(qq^gG) = F 1 {t/s,u/s), (3) 

dt v ; 36sMp u ' ' 1 h K ' 

^(qg _> q G) = ^Q- F 2 (t/s,u/s), (4) 
dt v ; 96sM| w ' 1 h K ' 

^{gg _> gG) = F 3 {t/ S , u/s), (5) 

dt vyy y ; WsMp w 7 ; ' v ; 

where Q 2 stands for the characteristic energy scale of the hard scattering amplitude, q stands 
for the light quark u, d and s respectively, and the functions ^1,2,3 take the following form 

F 1 (x,y) = — (4xy + z)(l-2xy + z 2 ), (6) 
xy 

F 2 {x,y) = — {2x-z 2 -y 2 ){z + z 2 + x{4 + z)), (7) 

xz 

F 3 (x, y) = — [x 2 y 2 + 2xy{z 2 -z + l) + (l + z + z 2 ) 2 ] . (8) 

The relation z = 1 + x + y is implicitly adopted that can be deduced from the fact of 
s + t + u = m 2 . It can be easily find that ^1,3(0;, y) = Fx >3 (y, x) as is the requirement of the 
invariance under exchange of the Mandelstam variables t and u. 
For the background subprocesses ij — > Z°k, we obtain 

where xw = sin 2 ^, R q = — 2e q xw and L g = 1/2 — 2e g a;vi/ with e g stands for the electric 
charge of q- quark. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In the following, we shall first discuss the uncertainties in estimating the graviton/Z° 
production. For the LO estimation, we shall concentrate our attention on the main uncer- 
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tainties, which include the PDFs for the light quarks, the choice of the factorization scale 
fiF and the typical energy scale Q 2 for the hard scattering amplitude and the number of the 
extra dimension for the graviton production. 

The number of extra dimension (5) can either be too small or too large. We are interested 
in the case in which 5 is not too large (say 5^6), under such condition the mass splitting 
Am is so small that the sum over the different Kaluza-Klein states can be replaced by a 
continuous integration, and then the enormous number of accessible Kaluza-Klein modes 
can rightly compensate the 1/Mp factor in the scattering amplitude. While an even larger 
number of extra dimensions shall lead to the mass splitting Am become comparable with 
the experimental energy resolution, hence only a smaller number of Kaluza-Klein modes 
can be produced, and then the total cross-section is negligible due to the unavoidable 1 / Mp 
suppression. On the other hand, it is argued that 5 can not be too small as is required by the 
latest torsion-balance experiment 16(. Constraints from cosmology also lead to non-trivial 
bounds on extra dimensions 13, ll8[ 1 . So, we take 5 £ [2,6] to do study its effect to the 
hadronic production. 

As is shown in Eq.(TjQ), PDFs F l Hl Pl (xi, fi 2 F ) and Ffj 2 p 2 (x2, H%) generate certain uncertain- 
ties in the estimation. PDFs are of non-perturbative nature, and in Eq.([T]) they have been 
factorized out at the energy scale fMp with the help of pQCD factorization theorem. The 
PDFs can be obtained only through global fitting of the experimental data and evolute them 
to the requested characteristic scale in standard way of pQCD, so there are several groups, 



e.g. CTEQ [211], GRV [22j and MRS [23[ etc, who devote themselves to offer accurate PDFs 
to the world, and to keep PDFs updated with the newly available relevant experimental 
data. Thus in literature, different versions of PDFs (including different issues by the same 
group) are used in the estimates of the hadronic production. To be self-consistent with the 



LO pQCD calculation, we shall adopt the two LO PDFs: CTEQ6L [2jJ and MRST2001LJ23] 
as typical examples for PDFs. The versions of the gluon distributions ended with 'L' are 
accurate up-to the leading logarithm order (LLO), i.e., their QCD evolution effects with a s 
running are included, so for the production to show the uncertainties correctly up to LO 
accuracy, it is necessary for the PDFs, the hard subprocess and the QCD 'coupling constant' 



1 By using a hyperbolic curved other than the flat extra dimensions |19j, it has been argued that those 
cosmology constraints can be naturally satisfied Q. 
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a s 'run' to the energy scale Q 2 properly. When computing the production and taking the 
PDFs from one version of the three groups, the running a s should also be taken from the 
same group. 

As for the leading order estimation, how to choose the energy scale Q 2 is a very tricky 
problem. If Q 2 is chosen properly, the results may be quite accurate. From experience, 
for a hard subprocess with two-body final state, generally the choice of Q 2 = s/4 can lead 
to an accurate LO result. To see the uncertainty caused by different choices of Q 2 , we 
take two typical types for Q 2 , i.e. Type A: Q 2 = s/4 with s the squared center of mass 
energy of the subprocess; Type B: Q 2 = M 2 = p 2 + m 2 , the squared transverse mass of the 
graviton/Z° respectively. And further more, To see the uncertainties from Q 2 choice, instead 
of variation on the choices with Q 2 = /if,, the authors in literature, such as Ref. 24j, also try 
Q 2 7^ /if and see the uncertainty. Here following them, we also calculate the distributions 
with Q 2 ^ /if,. More explicitly, as suggested in Refs. [24 l25j . we take /if G [M 2 /A, AM 2 ] and 
Q 2 e [M 2 /A, AM 2 } simultaneously to do the discussion. 

Other parameters like m^o, x\y and the fraction of Z° decaying to invisible particles 
affect the production slightly, so we directly take them to be their center values as adopted 
in the literature, e.g. m z o = 91.187GeV, x w to be 0.2311 and the fraction of Z° decaying to 
invisible particles like (anti) neutrinos v{y) to be the center value of (20.00 ± 0.06)% 15]. 

Another important thing we need to be careful is the regulation of the cross section to 
avoid the non-perturbative regime, i.e. to deal with the signatures arising from collisions 
with parton center-of-mass energies of order Mb or above properly. At Vs ^> Mjj, parton 
collisions are expected to produce classical black holes 26J]. As has been pointed out that 
the discrepancy between the two cases with or without the cut s < M D shall be increased 



with the increment of E™ l -l t and decreased with the increment of Mo 12], therefore by 
taking proper Erjet cut in experimental analysis, one can select different ranges of Me, to 
probe perturbatively. And numerically, we can find that under the present conditions, with 
or without the cut s < M D shall affect our final conclusions slightly, so we shall adopt the 
cut s < M D in all the following discussions. 
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FIG. 1: Total jet+nothing cross-section at LHC integrated for the requirement that £?r,jet > E*t Jet 
with an acceptance cut on the jet rapidity \rj\ < 3. The dash-dotted line is SM background, the 
solid, dashed and dotted lines are for 5 = 2, 4 and 6 respectively. 

A. The uncertainties in estimates 

Firstly, we show the uncertainty caused by the different choice of extra dimensions 5. 
For the purpose, we fix PDF to be CTEQ6L, M D = 5TeV, Q 2 = s/4 and n% = Q 2 . We 
show the signal and the background rates for the transverse jet energy larger than E™™ t in 
FIGJH with an acceptance cut on the jet rapidity \rj\ < 3. It is shown that the signal rates 
decrease with the increment of 5. The SM background is bigger than the graviton signal in 
the lower transverse jet energy region but it drops down much quickly than the signal, and 
one may distinguish the signal from the background in large transverse jet energy region. 
In another words, the large transverse jet energy region shall provide an effective platform 
to distinguish the signal and the background. 

Secondly, we show the uncertainty from different choices of Q 2 and /if,. For such purpose, 
we fix PDF to be CTEQ6L, Mo = 5TeV, 5 = 4. It is found that by taking two choices 
of Q 2 (type A and type B) under the case of fi 2 F = Q 2 , the uncertainties for both the 
background and the signal are small, i.e. the differences between these two types of Q 2 are 
less than 10% for both the background and the signal. Further more, we show the case of 
(Ap 7^ Q 2 by shaded band in FIGJ2J with an acceptance cut on the jet rapidity \rj\ < 3, 
where /i 2 F e [M?/4,4M?] and Q 2 G [M t 2 /4, 4M t 2 ]. It is found the largest value is obtained 
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FIG. 2: Total jet+nothing cross-section at LHC integrated for the requirement that E^jet > E F 
with an acceptance cut on the jet rapidity \rj\ < 3. The thinner and thicker shaded bands stand 
for the SM background and the signal respectively, with the upper edge for \j? F = M 2 /A and 
Q 2 = M 4 2 /4, and lower edge for \j? F = AM? and Q 2 = AM?. 



TABLE I: Total jet+nothing cross-section a(pp — > jet + J3t jet) ( m unit: fb) versus PDFs at LHC 
integrated for the requirement that Exjet > E™ 1 ^ with an acceptance cut on the jet rapidity 
\r]\ < 3. Two LO PDFs: CETQ6L and MRST2001L, and one NLO PDF: CTEQ6M, are adopted. 





Signal with E T ,jet > E^f et 


SM background with E TJet > Eff et 


PDFs 


0.50 TeV 


1.0 TeV 


1.5 TeV 


0.50 TeV 


1.0 TeV 


1.5 TeV 


CTEQ6L 


32. 


5.1 


0.92 


46. 


1.1 


0.070 


MRST2001L 


31. 


4.5 


0.73 


47. 


1.0 


0.060 


CTEQ6M 


40. 


6.3 


1.1 


57. 


1.1 


0.074 



when fx 2 F = M 2 /4 and Q 2 = M 2 /4, and the lowest value is obtained when \j? F = AM 2 
and Q 2 = AM 2 . And from FIG [5] one may also observe that while the uncertainty for 
the background changes to be around 20% (as shown by the thinner shaded band), the 
uncertainty for the signal can be changed up to 50% (as shown by the thicker shaded band). 
Thirdly, we show the uncertainty from different choices of PDFs by fixing Mp = 5TeV, 



S 



100000 
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■ 
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FIG. 3: Total jet+nothing cross-section at LHC integrated for the requirement that Exjet > -^r jet 
with an acceptance cut on the jet rapidity \rj\ < 3. Two typical LO PDFs, e.g. CTEQ6L, and 
MRST2001L are adopted. 

5 = 4 and \j? F = Q 2 = s/4. We show the results for two LO PDFs in FIGJ3] with an 
acceptance cut on the jet rapidity < 3, i.e. CTEQ6L and MRST2001L. It is found that 
the results of CTEQ6L and MRST200L are close to each other, i.e the difference is less 
than 20%. More explicitly, we show the total jet+nothing cross-section a{pp — ► jet + fir, jet) 
versus PDFs at LHC in TABJU TABjT] shows more explicitly that even though the SM 
background is bigger than the graviton signal in the lower transverse jet energy region but 
it drops down quickly, and at E F l ^ et ~ ITeV it is less than 25% of the signal. At the present 
the next-to-leading order (NLO) results are not available 2 , and to have a rough idea on 
how NLO calculation will affect the present results, we take the 'miss matched' NLO PDF 



as CTEQ6M [21] to do the calculation. And it is found that the results from the 'miss 
marched' CTEQ6M shall be bigger than that of CTEQ6L by about 25%. So a full NLO 
estimation shall be helpful to improve our present estimations. 

Finally, we show the uncertainty from different choices of Mjj by fixing PDF to be 
CTEQ6L, 5 = 4, Q 2 = s/4 and \j? F = Q 2 . We show the signal and the background rates 
for the transverse jet energy larger than E F % ^et m FIGJH with an acceptance cut on the jet 
rapidity \{]\ < 3. It is shown that the rate decreases with the increment of Mb under the 



Such a NLO calculation is in progress [2 7j. 
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FIG. 4: Total jet+nothing cross-section at LHC integrated for the requirement that Ex jet > EJp J£ t 
with an acceptance cut on the jet rapidity \rj\ < 3. Three typical choices for = 3GeV, 5GeV 
and 7GeV are adopted. 

TABLE II: Corresponding Mp sensitivity ranges versus 5 with the high integrated luminosity 
C = 100/6 -1 . Mrj(3a) stands for the 3a exclusion limit and Mrj(5a) stands for 5a observable 
limit, where the center values are obtained by taking CTEQ6L, Q 2 = s/4 and \j? f = Q 2 and the 
errors are caused by the above mentioned main uncertainty sources. 



5 


Mi)(3a — exclusion) 


Mi)(5a — observable) 


2 


9.2 ± 0.7 TeV 


8.8 ± 0.9 TeV 


4 


6.1 ± 0.5 TeV 


5.9 ± 0.5 TeV 


6 


5.3 ± 0.3 TeV 


5.1 ± 0.3 TeV 



same condition Ex jet > Erjet- 

B. To what energy scale can LHC probe ? 

It is found that the SM background changes slightly within the reasonable regions of 
the above mentioned uncertain sources, so it can be treated as a basis to decide to what 
energy scale can LHC probe. Since the large extra dimensions can be probed only when the 
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deviation of the cross section within the framework of large extra dimension model from the 
SM background is large enough, we adopt the 5cr large extra dimension effect observable 



limit and 3a exclusion limit as suggested in the literature 14| to extract the constraint of 
fundamental energy scale Mjj, i.e. 



Act = a LE D ~ OSgd > £ (11) 

Aa = cr LED - a Bgd < . [L2) 

The corresponding sensitivity ranges versus 5 are shown in TABjTT]with high integrated 
luminosity £ = 100/6 -1 , where the cuts \r)\ < 3 and Exjet > 1-0 TeV are adopted. It 
has found that with integrated luminosity L ~ several fb~ l , the main uncertainty comes 



from the instrumental background 



0, 0, H, E 



9J, which includes both the systematic and the 



statistical errors. So we have taken a higher integrated luminosity C = 100/6 -1 to do our 
calculation such that the systematic error is dominant. And in doing the calculation, we 
require that a led > o'Bkgd, since it has been found that the present adopted effective gravity 



theory is mostly reliable in this region as shown in Ref.[12j. It is found that the effective 
scale M D can be probed up to (8.8 ± 0.9) TeV for 5 = 2, (5.9 ± 0.5) TeV for 5 = 4 and 
(5.1±0.3) TeV for 5 = 6, where the center values are obtained by taking CTEQ6L, Q 2 = s/4 
and ix 2 F = Q 2 and the errors are caused by the above mentioned main uncertainty sources 
that varies within their reasonable regions accordingly. 

Further more, one may observe that the center values for Mr, decreases with the increment 
of 5. Our present results for Md versus 5 with 5<r observable limit as shown in TABjTT] are 
consistent with the values of maximum determined in Ref.[12| (TAB. 3 there for the 
same integrated luminosity £ = 100/6 -1 ) within reasonable uncertainties, with the center 



value of our present one slightly bigger than that of Ref. 12], which is mainly caused by 



the fact that different criterion was adopted in Ref. 
about 10% is adopted to do the discussion. 



121 ]. i.e. a fixed systematic error that is 



It has been argued that 12| if the discrepancy for the results with or without the cut 
s < Mp becomes larger, then the ultraviolet contributions become important, and then our 
present estimation may be not under control. Numerically, we find that such discrepancy 
is small for 5 < 4 (e.g. it is less than 1% for 5 = 2 and 10% for 5 = 4) by taking the Md 
values listed in TABjTTl while for even larger 5 such discrepancy becomes quite large, e.g. for 
5 = 6 such discrepancy is up to 100% by taking M D = 5.1 TeV. This shows that our present 
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TABLE III: Maximum Md sensitivity versus S with the low integrated luminosity C = Wfb . 
The center values are obtained by taking CTEQ6L, Q 2 = s/4 and n 2 F = Q 2 and the errors are 
caused by the above mentioned main uncertainty sources. 



5 


2 


4 


6 


Max M D 


8.2 ± 0.5 TeV 


5.7 ± 0.3 TeV 


4.9 ± 0.2 TeV 



adopted effective theory may not be fully applicable for such a large extra dimension. One 
may hope to decrease such discrepancy by lowering the value of E™ l Jl t since a smaller E™ i p et 
leads to a smaller discrepancy, however by doing this, the probed Md sensitivity range shall 
only be slightly lowered due to the fact that E™ l ™ t can not be set too small otherwise it will 
be more difficult to distinguish the signal from the background, as has been shown in the 
last subsections that the background shall increase much more quickly than the signal with 
a decreasing E™jl t . 

It may be also interesting to make a discussion on the maximum Md sensitivity with 
smaller integrated luminosity, e.g. 10/fe -1 . Under such case both symmetric and statistical 
errors are comparable, and now the criterions ( TITIfTSl) are not applicable, so we adopt the 
criterion suggested by Ref. [12] to do the discussion, i.e. we add the two errors in quadrature 
and require 



a LED > V2-^ . (13) 

The corresponding Md sensitivity ranges versus 5 are shown in TAB IIIII with lower integrated 
luminosity C = lOfb^ 1 , where the cuts \r]\ < 3 and Erjet > 1-0 TeV are adopted and the 
errors are caused by the above mentioned main uncertainty sources that varies within their 
reasonable regions accordingly. 



IV. SUMMARY 



It is found that with sufficient luminosity at LHC the fundamental gravity scale can 
be probed up to several TeV, with the precise value depending on the number of extra 
dimensions. In the present paper, we have presented a detailed discussion on the leading 
experimental signal of real graviton at LHC based on the process pp — > jet +$t with the 
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help of the quantum gravity theory with large extra dimensions. The main standard model 
background to these processes together with their uncertainties have also been discussed. 
It is found that in higher transverse jet energy region, e.g. Erjet > 1-0 TeV, and with 
certain number of extra dimensions, the quantum gravity signal can be distinguished from 
the background. 
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